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A DIGITAL-TO- ANALOG CONVERSION CIRCUIT USING 


THIRD-ORDER POLYNOMIAL INTERPOLATION 

By William P. Dotson and Joe H. Wilson 
Manned Spacecraft Center 

SUMMARY 

This report describes the basic design for a circuit that is capable of reconstruct- 
ing digitized analog signals by interpolating between given sample points with a third- 
order polynomial. From an analytical description of the circuit inputs and outputs, 
design equations are developed and implemented on an analog computer. 

The theoretical error performance of the third- order digital-to-analog converter 
is developed from the analytical description, and the error performances are compared 
with experimental results. The same is done for a zero-order (sample and hold) 
digital-to-analog converter, and these results are compared with results of the third- 
order digital-to-analog converter. 

Both theoretical and experimental comparisons of third-order and zero-order 
digital-to-analog converters show a markedly superior performance of the third-order 
digital-to-analog converter. From these results, it is shown that the use of third- 
order digital-to-analog converters for analog-data reconstruction will allow a signifi- 
cant decrease in the bandwidth required for digital telemetry systems that transmit 
analog data, with no decrease in the quality of the data. 


INTRODUCTION 


In most telemetry systems in use today, digital-to-analog conversion is performed 
by zero-order (sample and hold) interpolation circuits. Such circuits have the attract- 
ive feature of relative simplicity. However, because of this attractiveness, the basic 
disadvantage of these circuits is frequently overlooked. This disadvantage is that 
higher than necessary sampling rates must be used in the digital encoding of the analog 
signal. The use of higher sample rates increases the bandwidth necessary to transmit 
the signal and thereby imposes additional cost and complexity on the remainder of the 
telemetry system. The design and performance of a digital-to-analog converter (DAC) 
that uses a third-order polynomial interpolation between sample points are presented 
in subsequent sections of this report; use of such a DAC permits lower sampling rates 
without decreasing the accuracy of the reconstructed analog signal. However, if the 
need to reduce the sampling rate (bandwidth) does not exist, the use of such a DAC will 
significantly improve the data quality. 



The theoretical and practical design and performance of a high- order (third order) 
interpolation circuit for digital- to- analog conversion of telemetry data are discussed in 
this report. The study described was initiated in an attempt to relieve the bandwidth 
problems in the NASA Manned Space Flight Network telemetry system that occur as a 
result of the high- sampling- rate biomedical data (i. e. , electrocardiograms (EKG's)), 
and, at the same time, to maintain the present accuracy level of the data. 


SYMBOLS 


A 

B, C, D, E 


e 

e in 

e out 
F, G, H 
f 
f 

c 

f 

s 

|G(f)| 

g(t) 

II, 12, 13 

K - V z i 

L 

m 


amplitude 

constant coefficients of polynomial terms 
feedback capacitance 

initial- condition voltage applied to the circuit 

voltage 
input voltage 

output voltage 

scale factors 
frequency, Hz 
cut-off frequency, Hz 

sampling frequency, samples/sec 

root- mean- square spectrum of g(t) 

time function 

integrators 


integrator gain, 1/ R-C^ 

ratio of the spectrum cut-off rate of the signal to 6 dB/octave 
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P', Q', R', S', T' waveform segments of an EKG signal 


P 

R 

R, 


R. 

1 


S 1’ S 2 


s 

T 

t 

TP1, TP2 

Y 

Y 

Y 

Y 

y(t) 

y'(t) 

z* 


p 

p 

T 

CO 


normalized sampling rate, samples/cycle 

resistance 

feedback resistance 

input resistance 

switches 

Laplacian operator 
period, l/f, sec 
time, sec 
test points 

reconstructed signal, V 

first derivative of Y with respect to time 

second derivative of Y with respect to time 

third derivative of Y with respect to time 

analog signal as a function of time, V 

analog reconstruction of y(t) 

feedback impedance 

input impedance 

error, percent or pulse code modulation (PCM) counts 

instantaneous error measured near the base-line crossing, 
percent 

instantaneous error measured at the peaks, percent 

reset time, sec 
sampling period, sec 
angular frequency, rad/ sec 



Subscripts: 

ic 


initial condition 


n — 0, 1} 2, 3, » . » 


THEORY 


The analog-to-digital encoding used in the Apollo telemetry system is described 
as follows and is illustrated in figure 1. 

1. The analog signal (fig. 1(a)) is sampled for its amplitude value at time inter- 
vals spaced r seconds apart (fig. 1(b)), where r is the inverse of the sampling rate. 

2. The amplitude values are then pulse code modulation (PCM) encoded 
(fig. 1(c)) before transmission. 

In zero-order digital-to-analog conversion, the PCM signal (fig. 1(c)) is accepted 
and converted to a pulse amplitude modulation (PAM) signal (fig. 1(b)), which is then 
fed to a holding circuit that provides an output curve similar to the curve shown in fig- 
ure 2(a). A third-order DAC will also accept PCM signals, convert them to PAM sig- 
nals, and then compute and output a third- order polynomial, which is a function of 
time, to fit the transmitted signal samples (fig. 2(b)). The reconstruction errors that 
result from this and other processes have been found analytically by previous work 
(ref. 1). Some of the results of that work, along with experimental checks performed 
in this study, are shown in appendix A. 

The purpose of this report is to describe a method of designing a circuit on the 
analog computer for the third- and zero-order DAC's described in reference 1. The 
geometries for calculating the errors for the third- and zero-order DAC's are pre- 
sented in figures 3 and 4. The use of the design method would allow an experimental 
verification of the error analysis shown in figures 3(c) and 4(c) and also would pro- 
vide the basic design for a prototype of the third- order DAC. 




(a) Analog signal. 


(b) Pulse-amplitude- modulation signal. 


Figure 1.- Analog-to-digital encoding. 
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(c) Pulse- code- modulation 
signal (eight quantiza- 
tion levels). 

Figure 1.- Concluded. 



(a) Geometry for calculating e^. 

Figure 3.- Reconstruction 



(a) Zero order. 



(b) Third order. 

Figure 2.- Zero- and third-order 
DAC outputs. 



(b) Geometry for calculating e . 

r 

errors for a zero- order DAC. 
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(a) Geometry for calculating e^. 



(c) Error as a function of normalized 
sampling rate. 

Figure 3.- Concluded. 


(b) Geometry for calculating e . 

Figure 4. - Reconstruction errors for a 
third- order DAC. 


The development of the design equations for a third- order DAC proceeds as fol- 
lows. The development is facilitated by initially assuming the existence of a PAM ver- 
sion of the analog signal, such as that shown in figure 5. At t = 0, the circuit must 
compute and output a signal that is described by a third-order polynomial that would 
connect the input PAM signal points , Yg , Yg , and Y^_ n - The output actu- 

ally starts at t = 0, n = 0 (or Yg) and continues only until t = r, n = 0 (or Yg), so 

that only the central two points of the given four are connected. This completes one cy- 
cle of operation. The next cycle begins at t = 0, n = 1 (or Yg) and continues until 

t = r, n = 1 (or Yg). This cycle could be visualized as sliding the data points leftward 

past a fixed time scale, so that time, for the circuit, is always between 0 and t. At the 
same time, the circuit, in an abstract sense, has both memory and precognition, be- 
cause the circuit must have knowledge of the samples Y^ and Y^ . 
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Error, c, percent 



(c) Error as a function of normalized 
sampling rate. 


Figure 4.- Concluded. 



Time.t 


Figure 5. - Third- order PAM version of 
the analog signal. 


This knowledge is necessary because 
the circuit output between t = 0 and t = r 
must be described by 


y(t) = Bt 3 + Ct 2 + Dt + E (1) 


which is a third-order polynomial. Because 
there are four unknowns (B, C, D, and E) 
in equation (1), the signal voltages (Yj, 

Yg, Yg, and Y^) must be given at four 

sample times (-r, 0, r, and 2t, respec- 

tively) (fig. 5). The coefficients B, C, 

D, and E may be evaluated from the follow- 
ing equation set. 


Y(-t) = B(-t) 3 + C(-t ) 2 + D(-t) + E = Yj 

( 2 ) 

Y(0) = B(0) 3 + C(0) 2 + D(0) + E = Yg (3) 


Y(t) = B(t) 3 + C(t) 2 + D(t) + E = Y 3 


(4) 


Y(2t) = B(2t) 3 + C(2t) 2 + D(2t) + E = Y 4 


(5) 
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As a result of solving equations (2) to (5) 


B = 


-Y 1 + 3Y 2 -3Y 3+ Y 4 
6t 3 


C = 


Y 1- 2Y 2 +Y 3 

2t 2 


D = 


- 2Y 1 - 3Y 2 + 6Y 3 - Y 4 
6t 


E = Y r 


The initial conditions of the output signal may also be calculated. If 

Y(t) = Bt 3 + Ct 2 + Dt + E 

Y(t) = 3Bt 2 + 2Ct + D 
Y(t) = 6Bt + 2C 
Y(t) = 6B 

then 

Y(0) = E 
Y(0) = D 
Y(0) = 2C 
Y(0) = 6B 


(6) 

(?) 

( 8 ) 

(9) 

(10) 

(ID 

( 12 ) 

(13) 

(14) 

(15) 

(16) 
(17) 
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I 


BASIC ANALOG DESIGN 


Equations (6) to (17) in the preceding section are sufficient for the design of an 
analog circuit that will perform the desired interpolation. The equations that specify 
the initial conditions (eqs. (14) to (17)) of the circuit are in terms of B, C, D, and 
E, which are in terms of Yj to Y 4 (eqs. (6) to (9)). The implication is that four 

samples, in proper time sequence, must be accessed simultaneously for computation 
of the initial conditions each time a new sample is received by the circuit. 

This delay can be readily implemented with a device similar to a shift register 

(fig. 6). The operation of this arrangement is as follows. If n = 0, then the registers 

hold, from top to bottom, Y„ = Y., Y„ =Y „, Y„ =Y„, and Y, =Y 1 . The next 

4-n 4 o-n o z-n z 1-n 1 

data value Y c will be read into the top 
5-n 



Figure 7.- Flow chart for the use of cir- Figure 8.- Flow chart for the use of cir- 
cuitry with digital logic elements. cuitry with analog logic elements. 
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The initial-condition circuitry, which is driven by the sample-delay circuitry, 
is formed of summing amplifiers designed according to equations (6) to (9) and 
amplitude- scaling constraints. Consideration of full-scale values for to and 

of various sampling rates indicates that the initial-condition circuitry will operate well 
if scaled according to the following equations, which are taken directly from equa- 
tions (6) to (9). The equations are developed in appendix B. 

The initial-condition circuitry should have the following four outputs. 

-T^=-I6(- Y l + 3Y 2- 3Y 3 +Y 4) < 18) 

-T = -S( Y 1 - 2Y 2 - Y s) <19) 

T=2?(- 2Y 1- 3Y 2 + 6Y 3- Y 4) (20) 

-E = -Y 2 (21) 


The initial-condition-circuitry design is 
shown in figure 9. (A more extensive dis- 
cussion of the use of the initial-condition 
circuitry can be found in appendixes C 
and D. ) Different scaling factors may be 
chosen, of course, provided appropriate 
rescaling in the following portions of the 
circuit is also completed. 

The remainder of the interpolation 
circuit consists of integrators that are 
driven by the initial-condition circuitry. If 
equations (17) and (18) are examined simul- 
taneously, it is clear that 


3 3 

- = ~ jq Y(0) (22) 




10 


Figure 9. - Initial-condition-circuitry 
block diagram. 



3 

Therefore, if -3 Bt /8 is input to an integrator with a gain of 
t^Y/ 8. This output requires an initial condition of 


-2/t, the output will be 


2 

T 

8 


Y(0) = 




(23) 


which is provided by the initial- condition circuitry, as specified by equations (16) and 
(19). 

2 2 

Equation (19) involves not +Cr , but -Ct , because integrators invert their 
initial condition. This procedure is continued through two more integrations, with 
initial conditions on each integrator initial- condition input appropriate to that integra- 
tor output. The final result is shown in figure 10. 



Figure 10. - Integrator -circuitry block 
diagram. 


The circuit discussed in this report 
requires three timing signals. One signal 
is a clock (or track) signal that is synchro- 
nous with the data. The other two signals 
(a hold signal and an integrator-reset signal) 
are derived from the clock signal. 

The track and hold signals are neces- 
sary only for the track and hold pairs that 
are used in a shift-register-like arrange- 
ment and may be eliminated if the sample- 
delay circuitry is implemented with digital 
logic. The track signal occurs at the leading 
edge of the clock pulse, and the duration of 
the track signal is 0. Olr. The integrator- 
reset signal must occur after the clock pulse 
so that the signals have completed shifting 
before the resetting of the integrators. The 
reset signal reinitializes the integrators 
each time a new sample value enters the 
circuit. The complete circuit is shown in 
figure 11 with track and hold amplifiers in 
the sample-delay circuit. Example analog- 
data and data-clock sequences are shown in. 
figure 12. 
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Figure 11. - Third- order DAC block diagram. 
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(a) Sample-and-hold reconstruction of (b) Data-clock (or track) signals 

an example analog- data word at TP2. 

sequence at TP1. 


Figure 12. - Example data and clock sequence. 







(c) Signal from figure 12(a) phased 
with the clock signal. 


(d) Signal from figure 12(c) delayed 
one sampling period. 
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Time, t 



(e) Signal from figure 12(c) delayed 
two sampling periods. 


(f) Signal from figure 12(c) delayed 
three sampling periods. 


Figure 12. - Concluded. 


EXPERIMENTAL RESULTS 


The digital-to-analog conversions made in this experiment were carried out on 
an analog computer. A function generator was used to obtain the sinusoidal and square- 
wave inputs, and a polyrhythm generator was used to obtain the complex EKG analog 
inputs. The outputs of the system were displayed on an eight-channel strip- chart 
recorder. 

To check the initial- condition- circuitry outputs, a square- wave input was used 
to simulate a step-function input. The initial -condition circuitry was then checked for 
zero outputs, because the mathematical derivative of a step function is zero. The 

functions Y, Y, and Y are shown to be zero in figures 13 and 14. The 0. 2-hertz 
input in figure 15 did not yield derivative values of zero because the sampling rate was 
only 5 samples/cycle, or 2. 5 samples/half-cycle. Four or more samples/half- cycle 
are required in order for the derivatives to settle to zero, because this would require 
that Yj = Y 2 = Yg = Y 4 , which implies that B = C = D = 0. This relationship can be 

seen by examining equations (6), (7), and (8). If B, C, and D are identically zero, 

then Y(t) = Y(t) = Y(t) = 0. (See eqs. (11), (12), and (13). ) These zero-derivative out- 
puts were obtained with no adjustments to the theoretical gains of the initial- condition 
circuitry. The integrator gains were checked with respect to the calibration proce- 
dures and were found to agree exactly with the theoretical values. 
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Figure 13. - Results of a 0. 05-hertz square-wave input 
(p = 20 samples/cycle) as a function of time. 



(c) Zero- order DAC output (1 sample/sec). I f I I I I ) ) I t I 4 



Figure 14. - Results of a 0. 1-hertz square-wave input 
(p = 10 samples/cycle) as a function of time. 
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Figure 15. - Results of a 0. 2-hertz square-wave input 
(p = 5 samples/cycle) as a function of time. 




To obtain the peak and base-line crossing errors, a comparison of the third-order 
DAC and the zero-order DAC was obtained by the method shown in figures 3 and 4. The 
results obtained with sinusoidal inputs of 0. 01, 0.02, 0.03, 0.05, 0.07, 0..1, 0.2, 0.3, 
and 0. 4 hertz are shown in figures 16 to 24. The error analysis is shown in figures 3(c) 
and 4(c). The accuracy of the results is ± 2. 5 percent. Because these data were ob- 
tained from the strip-chart recordings (figs. 16 to 24), the smallest scale obtainable 
was 5 volts. The experimental errors obtained were less than the theoretical errors 
in some cases, partly because of the ±2. 5-percent accuracy and the inability to re- 
produce exactly the base-line crossings and peaks of the theoretical analysis. As 
shown in figures 3(c) and 4(c), the theoretical and experimental errors closely agreed, 
and the general trends were the same. Because of processing delays within the cir- 
cuitry, the third-order DAC output lagged the input waveform by two sampling periods 
(2 seconds in this case) in figures 16 to 24. 

Figures 25 to 35 show a comparison of the EKG input from the polyrhythm gen- 
erator, the third-order DAC output at 100 samples/sec, the zero-order DAC output at 
100 samples/sec, and the zero-order DAC output at 200 samples/sec. The zero-order 
DAC output was artifically delayed two sample periods to permit accurate phase ad- 
justments for error calculations based on comparisons of both zero- and third-order 
DAC outputs with the original data. Because the input was from the polyrhythm gen- 
erator, the analog data were stable from cycle to cycle. Therefore, any variation 
in the DAC outputs was caused by reconstruction errors within the digital-to-analog 
conversions. 

Of the expanded versions of two of the inputs, the third-order DAC output most 
closely approximates the analog input (figs. 36 and 37). However, a slight error is 
present in the Q' wave in both figures on the third- order trace. The Q'-to-R' peak 
transition approximates a step function, and a step function with sharp corners requires 
an infinite number of terms to fit the curve. The third-order DAC generates only the 
first four terms of that series. The zero-order DAC is also in error at this point 
(fig. 37). 

It has been shown from both the analytical and the theoretical standpoints that a 
third-order DAC is more accurate than a zero-order DAC. The third-order DAC out- 
put is much smoother because the output is piecewise continuous. 

Figures 38 and 39 (cross plots of figs. 3(c) and 4(c)) show that, for normalized 
sampling frequencies above 5 samples/cycle, it is advantageous to use a third-order 
DAC instead of a zero-order DAC at twice the rate. This procedure results in a 
2: 1 advantage. This ratio becomes increasingly larger as p is increased for the 
zero-order DAC. This advantage is of major importance for two reasons. Use of the 
third-order DAC allows reduced sampling rates to be used with no degradation in ac- 
curacy and allows a substantial reduction in the bandwidth required to transmit the 
EKG's. 
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Figure 16. - Results of a 0. 01 -hertz sine-wave input 
(p = 100 samples/cycle) as a function of time. 






Volts Volts Volts Volts Volts ( Volts ( Volts 



Figure 17. - Results of a 0. 02-hertz sine- wave input 
(p = 50 sample^/ cycle), as a function of time. 


19 



Volts Volts Volts Volts Volts Volts Volts Volts 



Figure 18.- Results of a 0.03-hertz sine-wave input 
(p = 33. 33 sample$/cycle) as a function of time. 
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Volts Volts Volts Volts Volts Volts Volts Volts 



Figure 19. - Results of a 0. 05-hertz sine-wave input 
(p = 20 samples/cycle) as a function of time. 
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Figure 20. - Results of a 0. 07-hertz sine-wave input 
(p = 14. 3 samples/cycle) as a function of time. 



I 



Figure 21. - Results of a 0. 1 -hertz sine- wave input 
(p = 10 sample^/cycle) as a function of time. 
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Figure 23.- Results of a 0.3-hertz sine-wave input 
(p = 3. 33 samples/cycle) as a function of time. 
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Figure 25. - Sinus bradycardia. 



Figure 26. - Normal sinus 





Figure 27. - Sinus tachycardia. 



Figure 28.- First-degree block. 
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Figure 30. - Third-degree block. 



Figure 31. - Atrial flutter. 
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Figure 36. - Expanded view of normal sinus rhythm. 



Figure 37. - Expanded view of second- degree block. 
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figure 38. 


- Plot of sampling frequency required for a third- order DAC to equal the 
performance of a zero-order DAC at the base-line crossing. 



Normalized sampling frequency (zero-order DAC), p, samples/cycle 


Figure 39. - Plot of sampling frequency required for a third- order DAC to 
equal the performance of a zero-order DAC at the waveform peaks. 


CONCLUDING REMARKS 


The experimental results of sample data that were processed by the third- and 
zero-order digital-to-analog converters indicate a markedly superior performance on 
the part of the third-order digital-to-analog converter. The experimental results are 
in close agreement with the results predicted by the theory. 

It has been shown that the use of a third- order digital-to-analog converter will 
allow substantial reductions in the transmission bandwidth, with no decrease in the 
data quality. If no reductions in transmission bandwidth are necessary, then the use 
of a third- order digital-to-analog converter will greatly improve the quality of the re- 
constructed data. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, October 14, 1971 
921-10-00-00-72 
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APPENDIX A 
ERROR EQUATIONS 


The equations for the peak error e and the base-line crossing error e. for a 

P ^ 

zero- and a third-order DAC were derived in reference 1 and are shown as follows. 
The zero ordinate of each figure was shifted to make the error equations as simple as 
possible. 


ZERO-ORDER DAC 


Figure A-l is used to find the instantaneous error midway between two sample 
points for a zero-order DAC. This error is near the base-line crossing. The instan- 
taneous error at t = r/2 for the sample-point locations shown in figure A-l is 



Figure A-l. - Geometry for calculating 
for a zero-order DAC. 
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where p = normalized sample rate in samples/cycle 
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Amplitude 


(A2) 


e b = 100 sin ^ 


Figure A-2 is used to find the instantaneous error at the peaks for a zero-order 
DAC. The instantaneous error at t = 0 for the sample-point locations shown in fig- 
ure A-2 is 



Figure A-2. - Geometry for calculating 

e for a zero-order DAC. 
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By using trigonometric identities, the equation 


e p = 100 ( 1 - cos 


(A4) 


can be obtained. Equations (A2) and (A4) are plotted in figure 3(c). 


THIRD-ORDER DAC 


Figure A- 3 is used to find the instantaneous error midway between two sample 
points for a third-order DAC. This error is near the base-line crossing. The instan- 
taneous error at t = 0 for the sample-point locations shown in figure A- 3 is 



Figure A- 3. - Geometry for calculating 
for a third-order DAC. 
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To find E, solve the set 


y, (- t ) “ B (- t ) 3 - c (- ¥) 2 + D (- 1) + E = A sin ["£ - T )] (A7a) 

y '(- 5 ) - B (- 5 ) 3 + C (- if + d(- i) + E = A sin (» D - 0 < A7b > 



(A7f) 


Because <u = 2ir/T and r = T/p 



(A7g) 

By using trigonometric identities, the equation 

e, = 10o[- sin - + i flO sin — - sin — ) (A8) 

b |_ P 16 \ P P / 


can be obtained. 



Figure A- 4 is used to find the instantaneous error at the peaks for a third- order 
DAC. The instantaneous error at t = 0 for the sample-point locations shown in fig- 
ure A- 4 is 



Figure A-4. - Geometry for calculating 

e for a third- order DAC. 
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To find E, solve the set 

y, (- x) = B (' x) + C (“ x) + D (- x) + E - A Sin [a>d - (Alla) 

y, (-j) =B (-i) +c ("5) + D (- i) + E = A sin ["(x ~ (Allb) 

+c © + D S) + E = A sin [ w (J + ^)] (Aiic) 
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y’(r) ' Kt) + C (t) + D (¥) +E = Asln[o,(J+^| 


(Alld) 


which results in 

E = TG ( 9 { sin [ w (? “ 5 )] + sin [ w (t + 1)]| 




(Alle) 


where co = 2?r/T and t = T/p. By using trigonometric identities, the equation 


e = 

P 



7T 

cos 

p 



(A12) 


can be obtained. 

A comparison of the errors obtained with a third- order DAC and a zero- order 
DAC can be determined by cross plotting the normalized sampling rate of each and in- 
dicating the percent errors at various points, as is done in figures 38 and 39. The 
figures show that for normalized sampling frequencies above 5 samples/cycle, it be- 
comes advantageous to use a third- order DAC instead of a zero- order DAC. In fig- 
ure 39 — the worst case of the two — for p = 10 on the zero-order DAC, p « 5 on 
the third- order DAC for the same accuracy. 
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APPENDIX B 

SPECIFIC DESIGN FORMULAS 


The basic design element in analog circuitry is the high- gain amplifier, which is 
indicated by the following symbol. 




e out 

— o 


Ideally, the high-gain amplifier has infinite open- loop gain, infinite input impedance, 
and zero output impedance. Normally, the gain is also assumed to be negative; that 
is, the amplifier inverts the signal. Under these conditions, for the following circuit 
configuration, equation (Bl) is applicable. 



*° e out 


6 j. = - ttt e. = -Ke. 
out Z. in in 
1 


(Bl) 


The gain for this circuit configuration, then, is a function of only Z ^ and Z.. 

Rather than always drawing in the impedance elements, a new symbol is usually 
defined 



where K is the ratio of the feedback impedance to the input impedance. Summing am- 
plifiers are easily constructed from the previous configuration and have output voltages 
given by equation (B2). 



Integrators are formed by changing the feedback element to a capacitor and have 
output voltages given by equation (B3) . 



(B3) 


The term l/s is a Laplacian operator that indicates integration. The term l/R.C^ 

controls the gain involved in the integration. The previous configuration may be 
abbreviated to 


where L = l/R.C f . 
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Integrators can also be configured to form the integral of summed inputs. 



Initial conditions must be known before solutions to integrals can be found. Known 
initial conditions can be placed on integrator outputs by appropriately charging the 
capacitor in the feedback loop. The time in which this capacitor is being charged is 
called the reset time. Integrators, then, have additional circuitry that is used for 
setting or resetting the initial conditions on output. One possible circuit follows. 



Because of the assumed infinite gain of the amplifier, the point A is at a virtual 
ground; therefore 


^ic e out c d ( e out) _ 0 


(B4) 


Equation (B4) can be written as follows by using Laplace-transform notation. 


^ic e out 


sR. R. 

1C 1C 


+ C.fse - e ,(0)1 = 0 
f|_ out out' J 


(B5) 


44 



or 


'out 


fe +ac <h4 +c ^ (0) 


<W<°> A 


'out 1 

s + 


R ic C f 


^ic r i 

' R ic c fppn 


nS) 


e ogt (t > “ e out (0)e 


y R ic C f Ac 1 1 - a' t/Rli:Ct l 


R ic c iV 


e out (,) = e out (0 > e 


t/R. C, ( -t/R. C,\ 

/ 1C f - sAl -e / 1C V 


(B6) 


(B7) 


(B8) 


(B9) 


It should be apparent that a definite time interval is required to set the initial 
conditions on the integrator output. At t = 5R. c Cj, the output will have reached 

99-1/3 percent of its correct value. After this time has elapsed, the integrator may 
be switched to the operate mode. The reset time then should be 


p s 5R- c C f 


(BIO) 


For a duty cycle of 99 percent (operate- to- reset ratio), p^ 0. Olr, where t is the 
time interval between samples. A further constraint on the integrator circuitry is in 
its gain. The integrator gain L is 




(Bll) 
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APPENDIX C 

SCALING COMPUTATIONS 


The scaling equations are developed by conjecturing a normalized full-scale 
sine- wave input sampled at a very high rate so that, in the limit, the sequential values 
of E (fig. 11) approach Y. That is 

lim E *■ Y = sin cot (Cl) 

p — 0 ° 

The first three derivatives of Y can then be found. 

Y = co cos cot (C2) 

Y = -cousin cot (C3) 

Y = -oo^cos cot (C4) 

If equations (8) and (15) from the text are equated 

Y(0) = D = ^ (-2Y X - 3Y 2 + 6Yg - Y 4 ) - co cos cot (C5) 

If equation (C5) is manipulated and multiplied by the scale factor F 

F6tY( 0) = F6rD = F(-2Yj - 3Y g + 6Y 3 - Y 4 ) - F6rco cos cot (C6) 


If sin cot is a full-scale output, then F6rco is chosen such that F6rco cos cot is one- 
half full scale. (One-half full-scale output is chosen arbitrarily to avoid both of two 
common error sources in analog devices, which are (1) the nonlinearity errors at full- 
scale voltage outputs and (2) noise errors at voltage outputs near zero. ) Therefore 


F6tco = ^ 


(C7) 
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and 


F = - 

* 12 TO) 


(C8) 


Because o> = 277f, t = 1/f , and p = f /f, 

s s 


F = 


P 

2477 


(C9) 


Equation (C9) implies that the scale factor depends on the lowest samples/cycle de- 
sired. If a rate of only 77 samples/cycle is desired, p = 77 and F = 1/24. The 
initial condition on Y should then be scaled to be 


6tY(0) _ tY(0) _ Dt 
24 " 4 ~~T 


(CIO) 


so that 


T 1 = Yi (‘ 2Y 1 - 3Y 2 + ff 3 - Y 4) 


(Cll) 


Similarly, the scale factor G on Y(0) can be found. 


Y(0) = 2C=4(Y 1 -2Y 2+ Y 3 ) 


(C12) 


where 


lim -4 (Y. 

p — 00 T z ' 


2Y 2 + 


Y s) — 


-Cl) Sin Cl) t 


and 


(C13) 


Gt 2 Y(0) = 2Gt 2 C = G(Y X - 2Y 2 + Yg) 


(C14) 
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where 


lim 

p — 00 


°( Y 1 - 2Y 2 + Y 3) 


2 2 

- Gt co sin cot 


Because 


Gt 2 co 2 = 1/2 


G = 


o 2 2 

2co T 


and 


For p = n, G = 1/8. Thus 




Gt 2 Y(0) = 


t 2 Y(0) _ Ct 2 
8 "4 


and 



= l( Y i-* Y 2 + y a) 


Similarly, the scale factor H on Y(0) can be found. 

Y(0) = 6B = \ (-Yj + 3Y 2 - 3Y 3 + Y„) 

T 


where 


lim 

p — o° 



+ 3Y 2- 3Y 3 +Y 4) 


-CO 


3 


cos cot 


(C 15) 


(C16) 


(C17) 


(C18) 


(C19) 


(C20) 


(C21) 
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and 


Ht 3 Y(0) = 6Ht 3 B = H^-Yj + 3Y g - 3Yg + Yj 


where 


lim Hf-Y, + 3Y 0 

p — 00 ' 


3Y 0 + Y 


4) 


3 3 

-Hr co cos cot 


Because 


„ 3 3 
Hr co 


1 

2 


and 


H = 


3 

P 

16? 


for p = 77, then H = 1/16. Thus 


Hr 3 Y(0) = ^ Y(0) = 2gf- ’ T6 (" Y l + 3Y 2 ‘ 3Y 3 * Y 4 ) 


where 


3Bt 3 

8 



+ 3Y r 


3Y 3 +Y 4) 


(C22) 


(C23) 


(C24) 


(C25) 


(C26) 


(C27) 
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APPENDIX D 

CALIBRATION PROCEDURES 


The initial- condition circuitry shown in figure 9 computes numerical approxima- 
tions to successive derivatives of the original waveform. The derivatives are then 
used to drive the integrator circuitry. Obviously, if the derivative calculations are in 
error, the final output will also be in error. 

One means of calibrating the initial- condition circuitry is to put a sampled func- 
tion with known derivatives into the total circuit and then examine the outputs of the 
initial- condition circuitry. The easiest function to work with in this instance is a step 
input, so that in the steady state, all derivatives are zero. With this input, the outputs 
of the initial- condition circuitry should be identically zero; if they are not, the potenti- 
ometers in front of the summing amplifiers require adjustment. 

The integrator gains may also require adjustment. Their theoretically required 
gain 1/t may not be realized because of tolerance limits on the feedback capacitors, 
input resistors, et cetera. For this reason, potentiometers were included in front of 
each integrator for fine adjustments. 

From an inspection of the circuit shown in figure 11, it should be clear that the 
integrator chain serves only to fill in estimates of data values within the given sequence 
of data values Yg, which is the initial- condition input of the final integrator. It is in- 
tuitively obvious, then, that the potentiometers in front of each integrator should be 
adjusted so that the final output is piecewise continuous. A sampled sine- wave input 
is easiest to work with in this instance. 

The outputs of the three integrators can be checked to see if they are piecewise 
continuous as follows. First, the output of the integrator will be determined. The 
output of the first integrator (II) can be determined by 
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The output can be checked to determine if it is piecewise continuous by evaluating 

and comparing the result with the initial condition at the next set of data values. 

That is, the continuity can be checked by asking whether 


3Bt 2 Ct 2 _ C’t 2 
— 4— *+-4-- 4 


(D3) 


After solving, it is determined that 


Y 2- 2Y 3 +Y 4 = Y 2- 2Y 3 +Y 4 


(D4) 


Therefore, the output of II is piecewise continuous. 

The output of integrator two (12) can be found as follows. 


12 




2 Dr 
dt -T" 


(D5) 


or 


TO 3Bt .2 Ct . Dt 
I2 = - — t 


(D6) 


The output of 12 can be checked for piecewise continuity by asking whether 


12 


k-r ’ - 


D'r 


(D7) 


or 


3Br 2 Ct 2 Dt 


D’t 

4 


(D8) 
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When equation (D8) is solved, it is found that 


-Yj + 6Y 2 - 3Y 3 - 2Y 4 ft 2Y 2 + 3Y g - 6Y 4 + Y 5 (D9) 

Therefore, the output of 12 is not necessarily piecewise continuous. 

The output of integrator three (13) is 

13 = -/(• t t2 - ¥ 4 - ^)(?) dt + e (dio) 

or 

13 = Bt 3 + Ct 2 +Dt + E (Dll) 

The output of 13 can be checked for piecewise continuity by determining whether 

I3| t _ r = E’ (D12) 


or 



♦ (-«1 - « 2 + SY 3 - y 4 )(<f) + Y 2 = Y 3 


(D13) 


When equation (D13) is solved, it is found that Y g = Yg. Therefore, the output of 13 
is piecewise continuous. 
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